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Abstract 

The cortical magnification matrix M is introduced founded on a notion similar to that of the scalar cortical 
magnification factor AI. Unlike AI, this matrix is suitable to describe anisotropy in cortical magnification, which 
is of particular interest in the highly gyrified human cerebral cortex. The advantage of our tensor method over 
other surface-based 3D methods to explore cortical morphometry is that M expresses cortical quantities in 
the corresponding sensory space. It allows us to investigate the spatial relation between sensory function and 
anatomical structure. To this end, we consider the calcarine sulcus (CS) as an anatomical landmark for the 
primary visual cortex (VI). We found that a stereotypically formed 3D model of VI compared to a flat model 
explains an excess of cortical tissue for the representation of visual information coming from the horizon of the 
visual field. This suggests that the intrinsic geometry of this sulcus is adapted to encephalize a particular function 
along the horizon. Since visual functions are assumed to be M-scaled, cortical folding can serve as an anatomical 
basis for increased functionality on the horizon similar to a retinal specialization known as visual streak, which is 
found in animals with lower encephalization. Thus, the gain of surface area by cortical folding links anatomical 
structure to cortical function in a previously unrecognized way, which may guide sulci development. 



1 Introduction 

The patterns of the highly folded surface of the cere- 
bral cortex are prominent features of the human 
brain (Fig. 1 A). Primarily, folding permits a larger 
cerebral cortex surface area to fit inside the skull. 
However, folding ensures that the additional surface 
area is not homogeneously distributed, if the sur- 
face becomes intrinsically curved. The surface gain 
is spatially concentrated in certain cortical regions. 
The cortex represents sensory information in distin- 



guishable fields and therefore the question of the re- 
lationship between anatomical structure and sensory 
function is naturally given. We utilize methods of 
continuum mechanics and complex analysis to ex- 
plore this relationship in the cortex. 

Many studies of human cortical architecture 
show that sensory and motoric fields have some re- 
lationship to the gross sulcal and gyral morphology, 
although a substantial variability in both the size 
and location is observed [TIB]- In a few cases very 
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precise correlations between sulci and functional en- 
tities could be demonstrated. Motor cortex can be 
identified by the position of the central sulcus [5] and 
the primary auditory cortex has a clear spatial rela- 
tionship with Heschl's gyrus [HU?]. The most reliable 
relation is, however, the calcarine sulcus (CS) as a 
landmark of the primary visual cortex (VI) [SUlQ], 
Anatomical identification is also quite reliable for vi- 
sual areas outside VI, e. g., V5 lies at the intersection 
of the ascending limb of the inferior temporal sulcus 
and the lateral occipital sulcus [TTlll2j . 

In this paper, we consider VI because not only 
its structural anatomy but also its functional retino- 
topy, i.e., the spatial organization of the neuronal 
responses to visual stimuli (see below), is well stud- 
ied in this field [IMHl- The CS, where VI is lo- 
cated, begins near the occipital pole on the medial 
surface of a hemisphere. It continues towards the 
posterior end of the corpus callosum (Fig. 1 A). We 
will present evidence that the 3D form of the CS 
indicates a selective magnification of the horizon of 
the visual field that is neither accounted for in stan- 
dard retinotopic maps nor refiected in the density of 
retinal ganglion cells, although a modest increase of 
retinal cell density, a so called visual streak, can be 
found [19H2I]. 

Already in 1984 Rovamo and Virsu [22] have 
noted that an locally isotropic (independent of the 
direction of measurement) cortical magnification 
that is also symmetric with respect to the meridi- 
ans can be better approximated by taking the un- 
folded convex 3D form of the cortex into account. 
This does not imply that cortical magnification has 
to be strictly locally isotropic but curvature affects 
the overall layout. Recently, the influence of cor- 
tical folding in primate did also take into account 
the concave folds [53]. We apply tensor analysis to 
investigate how these symmetries, i.e., isotropy and 
meridional symmetry, relate to the gross folding pat- 
tern, in particular the concave fold of CS that cre- 
ates additional cortical space for the representation 
of the visual field close to the horizon. To this end, 
we compare intrinsically flat and curved surfaces of 
VI and investigate how the 3D form affects cortical 
magnification. We propose that in particular a hor- 
izontal stripe gains additional cortical space. The 
shape and location of this stripe is similar to the 
visual streak as a retinal specialization found to be 
very pronounced in some animals |24H30| . This sug- 
gests a link between cortical folding, M-scaling and 
the functional development of cerebral sulci. 



First evidence for the predicted selective cor- 
tical magnification of the visual horizon can be 
found in the literature: perceptual fiUing-in [31| and 
travelling migraine scotoma [32] . The big advan- 
tage of psychophysical methods that measure sco- 
toma is that neither method requires surface recon- 
struction to measure cortical magnification. But 
it is difficult to get reliable quantitative data, be- 
cause these methods involve psychophysical investi- 
gations with subjective evaluations from probands 
and patients. Furthermore, there is a brief report 
of such a phenomenon |33| . Only computationally- 
intensive methods, which allow precise surface-based 
morphometry using anatomical magnetic resonance 
imaging and functional magnetic resonance imaging 
(MRI and fMRI) [18l[3l[35], can provide a direct 
test of our predicted correlation between CS, as an 
anatomical cortical landmark, and increased M, as 
a functional cortical measure. 



2 Material and Methods 

We use concepts from continuum mechanics to de- 
scribe the retino-cortical map |36| . In the La- 
grangian description of continuum mechanics, the 
map is seen as a deformation rather than a transfor- 
mation. This description uses some coordinate sys- 
tem in the sensory space as reference configuration 
to express every quantity in the deformed configura- 
tion, i.e., as the cortical map of the sensory space. 
The deformation language is intuitive but needs to 
be adapted for our purposes. Wc emphasize that 
we take mostly a kinematic approach. No attention 
is paid to the dynamics of the creation of a retino- 
topic map other than that we compare two discrete 
stages of VI, one which is intrinsically flat and one 
that evolved from the latter such that it is intrinsi- 
cally curved with a major sulcus. The Lagrangian 
description is also suited for a dynamical description 
continuous in time of map development, for example 
with self-organizing neural networks |37| . 

Let us start by briefly clarifying the terminology, 
in particular related to cortical mapping and cur- 
vature; linear and areal cortical magniflcation; the 
naming of the visual coordinate system; and, finally, 
tensor versus matrix methods. 

Conformal maps allow us to describe cortical 
magnification by a scalar field, simply referred to 
as cortical magnification factor. In fact, we start 
to consider conformal maps between domains in the 
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Figure 1: (A) The calcarine sulcus (CS) begins near the occipital pole on the medial surface of a cortical 
hemisphere and continues toward the posterior end of the corpus callosum (thick black line). This is the 
site of VI. Its border is marked with a dashed line. (B) Medial side of the occipital pole, as in (A) but 
with straight course of CS. (C) The approximate 3D form of VI in the {u,v,l) coordinates (see text). The 
representations of the two vertical meridians are dashed and the fundus of CS is marked with a thick black 
line. (D) The retinal surface is approximated by a quarter of a sphere, or (E) by a flat disc, as assumed in 
many experimental data. (F) A Cartesian reference coordinate system of the visual field (^1,^2)- (G) and 
(H) Coronal section through VI. The red line traces the stria of Gennari. The M-L axis indicates the medial 
to lateral axis. 



3 



complex plane, that is, wc consider intrinsically flat 
domains. Conformal maps also exist between in- 
trinsically curved domains. Note that throughout 
the manuscript, we do not consider extrinsic curva- 
ture unless we say so. Of course, to fit larger cere- 
bral cortex surface area inside the skull, nature can 
make primarily use of extrinsic curvature, but not 
exclusively (unless the cortex would be a cylindri- 
cally scrolled structure). In the following, we will 
refer to VI as being curved ("curved VI"), if it is 
a surface with intrinsic curvature at least at some 
locations within VI, i.e., locations that have a non 
vanishing Gaussian curvature. As a consequence, a 
curved VI must be embedded in 3D even if it is 
mostly flat. And we will refer to VI as a "flat VI", 
if it is an intrinsically flat surface, i.e., the Gaussian 
curvature is zero everywhere. 

It can be important to distinguish linear and 
areal cortical magnification, if we assume Af-scaling, 
i. e., the fact that a measured quantity remains quali- 
tatively similar across the entire visual domain when 
magnified in inverse proportion to M . For example, 
the time required for a scotoma of certain size to fade 
and become replaced by its background (perceptual 
fiUing-in, see discussion) can be proportional either 
to linear or areal cortical magnification. The under- 
lying mechanism could involve merely the scotoma 
diameter, which would suggest linear M-scaling or 
involve the total area covered by the scotoma, which 
would suggest areal M-scaling. The term M-scaling 
is unfortunately sometimes used without mentioning 
which is meant. Note that the distinction between 
linear or areal cortical magnification is in cases where 
Al is a scalar field not essential as these factors can 
easily be converted, but one still must specify which 
method was used and in which units M is given (see 
next section). 

To specify loci in the visual hemifield, a polar 
coordinate system with its origin in the viewer's fix- 
ation point is used (Fig. 1 E). The coordinate lines 
at fixed polar angles are called meridians, with the 
horizontal meridian (HM) as the reference set to be 
zero degrees (bold solid line in Fig. 1 E). Meridians 
increase and decrease in the upper and lower visual 
field quadrant to 90 degrees going anticlockwise and 
clockwise, respectively, until the two vertical merid- 
ians (bold dashed lines in Fig. 1 E) bound the visual 
hemifield. 

We will express the generalized cortical magnifi- 
cation in terms of a matrix M . We want to stress 
that cortical magnification is a concept independent 



of the retinal and cortical coordinate systems chosen 
to represent it. Therefore, M actually is a cortical 
magnification tensor. It may lead to some confusion, 
when we later introduce M based on the cortical 
metric in the sensory coordinate system. We need 
to consider some sensory coordinate system, which 
merely refiects the fact that cortical magnification is 
not exclusively a property of the intrinsic geometry 
of the cortex but of the sensory map in the cortex. 
In mathematical terms, M is the square root of the 
matrix expressing the Riemannian metric of VI in 
the sensory polar coordinates. 

The Lagrangian formulation of a such maps is a 
tensor approach. As an alternative to Riemannian 
geometry, it provides by definition some retinal pa- 
rameterization of the cortical surface that is needed 
for the concept of magnification from retina to cor- 
tex. In contrast, the cortical metric in an arbitrary 
parameterization can be used, for example, to define 
the Laplace-Beltrami operator as a generalization of 
the Laplace operator and much effort has gone in 
the development of computationally advantageous 
parameterizations in this regard. However, since the 
cortical magnification is usually discussed referring 
to meridians and eccentricity, it is natural to use ex- 
clusively this polar coordinate system. Therefore, 
we refer to the generalized M as the cortical magni- 
fication matrix, and avoid the potentially daunting 
term tensor. 



3 Results 

We will first define the contour of an intrinsically fiat 
VI. Then we will define the cortical magnification 
matrix and use this concept firstly to compare cor- 
tical magnification along different radial directions 
(from horizontal to vertical) on the flat VI and show 
that the vertical direction has an increased corti- 
cal magnification. Let us emphasize that retinotopic 
maps that are assumed to be intrinsically flat can 
still be extrinsically curved. We will define an intrin- 
sically curved VI starting from an extrinsically fiat 
VI by plausible deformations. In a second step, the 
retinal coordinate system is defined on this curved 
VI, preserving some symmetry constraints while re- 
laxing others. Finally, we use the concept of the 
cortical magnification matrix to predicted a selective 
increase of cortical magnification along the horizon 
due to cortical folding. 
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Figure 2: (A) Retinal grid which serves as a reference frame in the Lagrangian formulation of retinotopy. In 
the inner 9°, iso-eccentricity lines are drawn in l°-steps in green. From 10° to 90° iso-eccentricity lines are 
drawn in 10°-stcps in blue. Meridians are drawn in 15°-stcps in red. (B) The neural representation of the 
retinal grid mapped with Eq. [3] serves as a flat model of VI. (C) The blue surface serves as a curved model 
of VI. It is created by forming B with lift functions (see Fig. 1 G and H) in the third dimension. The neural 
representation of the retinal grid on this surface is placed such that meridians are equidistantly spaced. 
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3.1 Symmetries in the retinotopic map 

As a starting point, wc need the contour of the area 
of VI. The contour of VI can be defined by the 
bounded domain of the visual hemifield that VI neu- 
raUy represents. Simply speaking, VI has a map of 
one visual hemifield with a nonuniform map scale. 
This scale is called cortical magnification. This map, 
or rather the retinotopic mapping function governs 
the contour of VI, because it maps the bounded do- 
main of the visual hemifield, approximated by a half 
disc, that extents to a visual eccentricity of nearly 
up to 90° (see Fig. 2A). We will contrast some em- 
pirical observations with deductive reasoning based 
on plausible symmetry conditions on the retinotopic 
map to obtain the contour of VI and to introduce 
these symmetries that guide us through the following 
sections. 

The retinotopic map is usually expressed by the 
cortical magnification factor Ad. For a mapping 
function from an one dimensional domain (a line) 
to another, e.g., along a single meridian, it is suffi- 
cient knowing the value of M on this domain. Note, 
that the term linear cortical magnification factor (see 
previous Section) refers to the one dimensional do- 
main. This can cause confusion, because M is often 
modeled as an inverse linear function of a single co- 
ordinate, the retinal eccentricity 0. A fact that will 
be important in the following brief analysis of the 
involved symmetries. 

In the one dimensional case M is simply the 
derivative of the mapping function. The mapping 
function can therefore be determined up to a con- 
stant of integration. This constant can be set to 
zero, because it describes only a translational shift. 

The complex logarithm provides a standard for- 
mulation of retinotopy in a fiat VI. Fischer 
first suggested this analytic function for the trans- 
formation of the visual field into its neuronal rep- 
resentation. While this relation was derived from 
visual inspection of ganglion cell density and recep- 
tive field size distribution assuming these quantities 
are M-scaled, it is supported by data of Ad [9) }39j . 
as shown by Schwartz |3D] with a power law \d~^ 
fit. The exponent n being sufficiently close to unity 
to be replaced by it. The meaning of parameter h 
will be explained below. Thus the integral equation 
reads as 

u{e) ^ j Ad{e)de = j ^Q-^de ^^\og9, (i) 

where u is a cortical linear distance, for instance the 



distance along the cortical representation of the hori- 
zon in the visual field. 

In two dimensions, the situation is more compli- 
cated. The exact unity value of the exponent n = 1 
can even be postulated ab initio from symmetry con- 
straints for 2D maps. It is therefore important to 
distinguish these two perspectives. In 2D, a possible 
but purely hypothetical function is obtained by gen- 
eralizing the real function in Eq.[l]to a corresponding 
complex function 

w = ^\ogz. (2) 

The magnitude of z is the retinal eccentricity and 
its argument <j> is the azimuth {z = 0e*"^), and the 
real and complex parts of w are Cartesian coordi- 
nates (u,w), respectively. To generalizing the real 
function in Eq. [T] to a corresponding complex func- 
tion implies two rather obvious constraints on the 
retinotopic map and one more subtle constraint. 

Firstly, a complex function implies conformal 
mapping, that is, cortical magnification becomes a 
scalar field. This is a symmetry, namely that cor- 
tical magnification at any location is invariant with 
respect to direction. Therefore, cortical magnifica- 
tion is locally isotropic. At least for some M-scaled 
visual functions this seems to be a natural symme- 
try requirement, such as for visual acuity. Secondly, 
complex functions are conformal maps between do- 
mains in the complex plane, that is, intrinsically fiat 
domains, which, of course, can still be extrinsically 
curved to accommodate the limited extend of the 
skull. This is another symmetry, namely that the 
Gaussian curvature is constant and zero everywhere. 

These symmetries, which are given by the con- 
straints of analytical functions, look rather reason- 
able. But even if these symmetries were reasonable 
in the light of missing or uncertain experimental data 
(see discussion), they should not be given preference 
to further symmetries that could also be found in the 
retinotopic map, in particular meridional symmetry 
of cortical magnification, i.e.. Ad is invariant unter 
retinal rotations around the center. 

We will discuss whether meridional symmetry of 
Af should be considered as a plausible constraint, 
but let us end this section by emphasizing its severe 
consequences. 

The real and imaginary parts of an analytic func- 
tion are conjugate harmonic functions that solve 
Laplace's equation. From Laplace's equation in po- 
lar coordinates, it can be easily shown that Eq. [2]is 
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harmonic and that any harmonic function that de- 
pends only on the eccentricity 6 must be of this form. 
Ahernatively, this follows from Cauchy-Riemann 
differential equations in complex analysis, which 
must be satisfied if we assume that the analytical 
mapping function is differentiable j41j . 

Thus, it is readily shown that any meridional 
symmetric analytic function must be of the form of 
Eq.O Therefore, any retinotopic map that is an an- 
alytic function with a meridional symmetric magni- 
fication factor Ad implies that cortical magnification 
is an inverse linear function (in fact, a linear func- 
tion with the axis intercept at zero). So if cortical 
magnification is investigated under the following two 
assumptions: (a) independent on direction (confor- 
mal map), (b) intrinsically flat (with (a) this leads 
to analytical functions), (c) meridional symmetry of 
M, one should be aware of the implicit assumption 
of inverse linear cortical magnification in the form of 
Eq.Hl 

Even more importantly, Eq. [5] fails to describe 
the retino-cortical map close to the representation 
of the fovea {9< 1) because the foveal point at 6' = 
is a singularity. To include the foveal region an offset 
can be introduced 

w = \\og{-z + l). (3) 
a 

together with another parameter a. Eq. [3] is not 
meridional symmetric. It follows from the reason- 
ing before that a retinotopic map that is based on 
an analytic function and includes the representation 
of the fixation point {9 = 0) cannot be meridional 
symmetric. 

Both parameters can now be interpreted easily 
when Eq. [3] is differentiated: ^ is the value of the 
linear cortical magnification factor M at the center 
of the visual field in terms of millimeters of cortex 
per degree of visual angle. The parameter b is the 
linear growth rate (slope) of the inverse linear cor- 
tical magnification factor on the real axis (0 = 0), 
i. e., along the horizon. This map is called monopole 
map (see Fig. 2 A-B). Values of the parameters a 
and b are given in the literature [12] • Note that 
sometimes different parameter are used, in particu- 
lar w = A\og{z + E2) with a ^ E2/A and b ^ 1/A. 

We use normalized values in units of a. The value 
of b is chosen as 0.57265a which would correspond to 
a = 0.117 and b = 0.067 as in Ref. [43]. When Eq.[2] 
is used, for example by |18| , there is only parameter 
b. Our dimensionless parameter a can be estimated 



from the contour of the flat map. 

The inverse linear magnification holds in the map 
given by Eq. |3]only at the horizon. Due to the shift, 
none of the other meridians take a simple course in 
the complex domain (overlapping with coordinate 
lines of the real and imaginary part or the coordinate 
lines of the absolute value and argument). There- 
fore, to investigate the cortical magnification factor 
along other meridians we introduce the cortical mag- 
nification matrix, which can express magnification at 
arbitrary points along arbitrary directions most eas- 
ily. Furthermore, unlike the scalar factor AI, this 
matrix is also suitable to describe anisotropy in cor- 
tical magnification, which is of particular interest in 
the highly gyrified human cerebral cortex. But it 
can also be useful for conformal maps, in which case 
the matrix at any point can be transformed to the 
identity matrix multiplied by a scalar. 

3.2 Generalized cortical magnification 

The natural description for a generalized cortical 
magnification is the Lagrangian description. It uses 
the retinal configuration as a reference to express 
every quantity in the cortical configuration. Mathe- 
matically this is expressed as 

$ : Dji Dvi : r 1-^ p = <l>(r), (4) 

where r and p are points in the retinal and cortical 
domain and Dyi, respectively. The Jacobian 
matrix of map $ can be interpreted as a deforma- 
tion gradient 

J<E,(r) = V$(r). (5) 

It is a homogeneous transformation tangent to the 
transformation $ attached at the image of r. Cau- 
tion is needed when the classical linear cortical mag- 
nification factors are derived from the components of 
J*. Usually non-Cartesian coordinates are utilized 
to describe the retinal reference state. Furthermore, 
contains a rotation which carries the retinal di- 
rections onto its neural representations. This rota- 
tion introduces problems if cortical magnification is 
not isotropic. To obtain the linear cortical magnifi- 
cation factor in accordance with its definition |39j . 
one must use the scalar product of a vector in the 
considered direction. The linear cortical magnifica- 
tion factor My^ of a retinal tangent unity vector Vr 
{\\vr\\ = 1) attached at r is the norm of its cortical 
image under the homogeneous transformation: 

M,^ = ||V$(«,)||. (6) 
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This norm and the underlying scalar product is de- 
fined via the cortical metric in retinal parameteriza- 
tion 

9vi = J^igRJiy (7) 

where gyi are the metrices of retina and cortex, 
respectively. The cortical metric Qyi leads to the 
definition of the cortical magnification matrix. This 
matrix is founded on a notion similar to that of the 
linear cortical magnification factor Af , but is much 
broader in conception. Once we have this matrix, 
we can calculate the linear and cortical areal magni- 
fication factors on the fiat and curved model of VI 
and compare the results. 

In continuum mechanics gyi is called the right 
Cauchy- Green tensor and termed C. For the sake of 
simplicity, we follow Salencon |36l p. 44] and adopt 
the term (cortical) expansion tensor for C. This 
name may better convey the meaning of this ma- 
trix. Its eigenvalues and eigenvectors give the shape 
and orientation of an ellipsoid representing an ini- 
tially spherical infinitesimal area in the retina. Fur- 
thermore this alternative name reminds us that we 
are not interested in the cortical metric expressed 
in an arbitrary coordinate system but in one that is 
related to the sensory space. This would be suffi- 
cient for other surface-based 3D geometric tasks us- 
ing, for instance, the Laplace-Beltrami operator to 
create flat maps with minimal surface area distor- 
tions. Using this particular cortical surface param- 
eterization, namely some sensory coordinates makes 
the choice fundamental. It links the metric to sen- 
sory function and it allows us to describe changes of 
the cortical surface from the sensory perspective. 

The square root of C is a pure stretch tensor that 
is directly related to the classical concept of cortical 
magnification and it could therefore also be termed 
cortical magnification tensor, yet the term matrix is 
probaly more intuitiv in the neuroscience community 

M = VC. (8) 

The matrix M is alternatively obtained from a po- 
lar decomposition of J^J^ji. The Jacobian matrix 
J^^ is the homogeneous transformation tangent to 
^f, which maps Cartesian coordinates to the visual 
field coordinates, i.e., either spherical or polar (reti- 
nal) coordinates. The distinction between spheri- 
cal and polar coordinates should be formally made, 
since data of the visual field position can be obtained 
either with perimetry or campimetry. We can avoid 
this, however, if we transform the retinal system into 



the Cartesian coordinates (^1,^2) of the visual field 
(see Fig. 1 G), because data is often available in 
these coordinates, for instance migraine aura sco- 
toma drawn on paper and subsequently scanned. 

3.3 Linear cortical magnification factor derived 
from the matrix M 

The linear cortical magnification factor M-d^ can be 
defined at any point r for arbitrary directions v. In 
experiments, M-d^ is usually measured along a con- 
stant meridian, though not necessarily along this di- 
rection. In conformal maps, e. g., the monopole map 
(Eq. [3]), depends only on position r but not on 

direction v. 

We compare the M^^, on the monopole map (flat 
VI) with the one on the curved VI along the 9- 
direction. This direction is most frequently used in 
the definition of the classical linear cortical magnifi- 
cation factor. Therefore, we will use M without any 
index for linear magnification in this direction. It is 
convenient to plot the inverse linear cortical magnifi- 
cation factor, because this is a linear function in 9 on 
HM in the monopole map. As the meridians in the 
monopole map change from horizontal to vertical, 
the inverse linear cortical magnification factor grows 
slower than linear (see Fig. 3). This is readily under- 
stood, considering the layout of the visual hemidisc 
mapped by Eq. [3] into its neuronal representation 
(Fig. 2 B). While the representations of all meridi- 
ans extend approximately equally long into the cor- 
tical u-direction (in anatomical terms the posterior- 
anterior axis), the more vertical they are, the larger 
is their evasion into the cortical f-direction (into the 
dorsal and ventral direction for the upper and lower 
visual quadrant, respectively). The increase in corti- 
cal space on the vertical meridian (VM) as compared 
to HM is about 20%. While in fact one study indi- 
cates the opposite asymmetry with HM being corti- 
cally over respresented with respect to the vertical 
meridian 1331. 



3.4 3D model of VI with curved retinotopy 

We construct the 3D model of VI based on our own 
and published data obtained at autopsy of neurolog- 
ically normal human brain [SlllOj. The dimensions 
of VI were determined from cross-sections in both 
hemispheres. The surface of VI can easily be identi- 
fied postmortem by the stria of Gennari. This band 
of myelinated axons can be traced (see Fig. 1 G-H). 
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Figure 3: Inverse linear magnification factor as a function of eccentricity 6 along different meridians (hor- 
izontal meridian (HM) 0°: blue; 30°: cyan; 60°: green; vertical meridian (VM) 90°: red, see inset). Solid 
lines are on the flat model, dashed on the curved model of VI. M has an inverse linear dependence on 
retinal eccentricity 6 only along the neural representation of HM on the flat model of VI (solid blue line). 
While on the flat model of VI, M increases with azimuthal distance at constant eccentricity (solid lines). 
This meridional asymmetry is reversed on the curved model of VI (dashed lines). In particular meridians 
close to HM gain cortical surface on the 3D model of VI. 
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Such traced curves provide lift functions l{u,v) for 
forming the originally flat monopole map (Eq. 
into its 3D from. We will only summarize the data 
obtained by visual inspection of the stria of Gennari 
used to define plausible functions that define the 3D 
form. In particular, we outline the restrictions we 
impose regarding plausible symmetry principles. 

First, we transform the complex function in Eq. 
[3] into the real domain 

//) ,\ 1 , /'a^+2b9cosU)a+b^9^\ 

u{e,4,) = 2fciog( ^ j, 

= iarctan(^j|§|;igL) . 

Eq. [5] describes together with the lift function /(m, v) 
a (0, (/)) -parameterized surface in M'^, that is, in the 
Cartesian coordinates {u,v,l). Note that the lift 
functions are used in a first step to describe the sur- 
face but not the retinotopy on it. We define lift 
functions in such a manner that the initially flat sur- 
face becomes intrinsically curved. Note that we may 
introduce a bias since we actually treat the intrin- 
sically flat surface also as extrinsically flat at this 
step of the construction. But the bias does not nec- 
essarily affect retinotopy because we rearrange the 
location of the meridians. Applying the lift func- 
tion directly to the retinotopic grid (Fig. 2 B) would 
drive some adjacent meridians farther away from 
each other than others, depending on There- 
fore, to obtain a curved retinotopy we need to re- 
arrange the location of the meridians by a inverse 
sampling technique, as will be described at the end 
of this section. 

VI is located entirely or nearly entirely on the 
medial surface of the occipital lobe. We align par- 
allel to this surface the {u, u)-plane and deflne the 
lift functions l{u,v) such that about two-thirds will 
lie within the CS walls. Futhermore, l{u, v) depends 
on the steepness of the walls of CS and on its course 
in the (u,w)-plane (see Fig. 1 A). Due to the large 
variations among individuals, we have to make some 
simpliflcations. We assume that the course of CS 
can be deformed to follow a straight course. 

To flnally construct the lift functions /(w, u), it 
is useful to deflne another landmark: the fundus 
of CS (FCS). Roughly speaking, it is the curve of 
maximum depth that spans the length of the CS. A 
more precise deflnition can be given based on cur- 
vature or based on distance functions [H]. Based 
on FCS, the assumption that the course of CS is 
straight can be formulated in a different way. In 
this case, the FCS assumes a curved line without 



torsion and its projection onto the (u, '(;)-plane is a 
straight line (see Fig. 1 C). For simplicity, we assume 
that CS has a bilateral symmetry with the symme- 
try plane going through this line. Furthermore, the 
mathematical description of VI is w.l.o.g. simplifled, 
if the projection of FCS onto the (u, ?;)-plane lies on 
the u-axis. Then the family of parametric profiles 
lu{v) = ^(u = const., v) in the coronal planes («,/)« 
completely describe the 3D-conflguration of VI. 

It significantly simplifies the later performed re- 
sampling to define the lift functions along curved 
coordinate lines in the (u, w)-plane with constant 6, 
instead of along the straight coordinate lines with 
constant v. Since the family of parametric profiles 
le {4>) live in curved planes we transform them into a 
Cartesian (^1,^2) plane (see Fig. 1 F). This results 

in a family of parametric profiles l^-i{S,2), which we 
finally choose to be Gaussian-shaped: 

hA^2)^de'^. (10) 

The parameter ct^^ defines the width of CS. Its value 
is (see Fig. 2 C). The parameter d gives the depth 
of FCS as a function of the new eccentricity coor- 
dinate ^1 in the {£,2,l)-p^s.nc. Again we choose a 
Gaussian-shaped profile 

= ^ |^l_e"%^ . (11) 

The parameter cr^j defines the steepness of CS. 
Its value is ^ (see Fig. 2 C). The parameter 
dmax = 170.94 gives the maximal depth of CS in 
units of a, equivalent to 20mm. The factor s ~ 

i J^^ \/ 1 + {%)^d6 is needed to resize the profiles, 
since we define the lift function on the (C17 C2) plane, 
but want to have a Gaussian-shaped FCS in the (u, I) 
plane. Each arc with constant eccentricity spans an 
angle tt in the visual field. On the (^1 , ^2) plane these 
arcs have all the same length, while on the (m, v)- 
plane their arc length is tts. Therefore, we have to 
resize ^^1(^2) by this quotient. 

The Gaussian curvature K is determined by this 
constrution because it depends only on the metric 
(the retinotopy depends on the metric in a retinal 
coordinate systems). In Fig. 4 we show isT in a color 
code on this surface. It is nearly everywhere close 
to zero, except for two main locations. One is at 
the saddle-shaped entrance of the FC. Here, K is 
negative (blue). Father along the fundus, when CS 
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reaches its maximal depth K changes to positive val- 
ues (red). 

To summarize, four hypotheses have been make 
about the area VI: (i) two-thirds will lie within the 
CS walls, (ii) the walls are symmetric with respect to 
the fundus, (iii) the shape of the walls can be approx- 
imated with smooth Gaussian-shaped profiles with 
constant standard deviation, and (iv) the shape of 
the fundus can also be approximated with smooth 
half Gaussian-shaped profile (we have alternatively 
chosen cosine-shaped profiles for both walls and fun- 
dus, with almost identical results) 



3.5 Retinotopy on a 3D model of Vl 

To compare retinotopy on the fiat and curved model 
of VI we now define the retinotopic grid on the lat- 
ter. Although this may seem obvious, we want to 
note, that we actually define the retinotopy on the 
curved 3D surface by this. 

Firstly, assume that HM maps to the fundus of 
the curved model surface. By definition HM splits 
the retinotopic map into two quarter fields, a dorsal 
{v > 0), lower field (</) < 0) and a ventral {v < 0), 
upper field {(j) > 0). If we simply lift the neural 
representations of the meridians from the flat sur- 
face model of VI onto the dorsal and ventral part 
of the curved model by Eq. [TOl the meridians drive 
away from each other in proportion to ^'g'^"^ . In 
fact, meridians should by equidistantly spaced in the 
curved model of VI, as they are on the flat surface 
(see red coordinate lines in Fig. 2 C). This can be 
achieved by rearranging (p on the flat model of VI 
before lifting with the inverse sampling equation 

0W = ^.-^(^^.(f)). (12) 

The function £e{(l>) gives the arc length of the 
profile le{(j)) between the origin (0 = 0) and the value 
of the argument (j). £g^ is the inverse function with 
the new value of (f) as the return value. The rear- 
ranged meridians are then lifted down from the fiat 
model with lift function (Eq. fTU)) . At last, we need to 
place the iso-eccentricity coordinates, i.e., the green 
(inner 10°) and blue (10°-90°) coordinate lines in 
Fig. 2 C. They are lifted from the fiat model with 
the lift function (Eq. llOp without rearrangement. 



3.6 Areal magnification and surface area gain 

Figure 5 A shows the areal cortical magnification fac- 
tor projected in the visual hemifield for the curved 
VI. This factor predicts the retinal ganglion cell 
density, assuming this quantity is M-scaled |20j . 
There is a larger value of det M on a stripe centered 
around HM. The inverse profile of M along this stripe 
is visualized in Fig. 3. The graph of M in Fig. 5 B 
visualizes this stripe as a shoulder in the exponen- 
tially decreasing curve AI{9) at HM. The shape and 
location of this stripe in the visual field (see Fig. 5 
A) is similar to a retinal specialization known as vi- 
sual streak. The visual streak is a stripe of elevated 
neuronal density along HM in the retina. It is found 
in some animals [24H30] but it is not pronounced in 
humans P"5H^ . The stripe centered around HM can 
be thought of as a virtual visual streak as described 
in the following. 

The neuronal density in the retina is assumed 
to be M-scaled (areal magnification). Let us define 
a virtual counterpart of the visual streak as the el- 
evation of inversely retinotopic projected densities 
of cortical cells. Along the horizontal meridian, the 
retinal representation of cortical cells is largely ele- 
vated with respect to other more vertical meridians 
outside this visual streak. A virtual visual streak is 
related to surface area magnification but also to sur- 
face gain by folding. Just as Eq.[6]defined the stretch 
of a unity vector, we can calculate the surface area 
gain of a unity area in the retina via the expansion 
tensor C. Considering rigid body transformations, 
C is the unity tensor 1 and surface area gain is a 
measure that must be zero. This led us to define 
the strain tensor £ = ^(C — 1). In continuum me- 
chanics s is called the Green-Lagrange strain tensor. 
The tensor e is a measure of surface gain from retina 
to cortex. The measure of surface area gain by cor- 
tical folding is obtained when the retinal reference 
frame is replaced by its neural representation on the 
flat model surface of VI. The strain tensor of the 
curved model of VI, with the cortical flat reference 
frame, is called E. It measures the surface gain by 
folding. 

The value of det E corresponds to the elevation 
of retinal density that is needed to satisfy M-scaling 
in a folded cortex. To show that this is in good agree- 
ment with some experimental data [HIIIM] (see also 
discussion), we look at det E along different merid- 
ians (Fig. 5 C). For all meridians, folding provides 
equally increased cortical space for the representa- 
tion of the periphery of the visual field. This in- 
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crease is mainly limited to the periphery because 
the calcarine sulcus just starts at the neural repre- 
sentation of the fovea and additional surface area is 
mainly gained when the calcarine sulcus has reached 
its maximum depth. Close to the fovea, meridians 
close to HM also gain cortical surface area for their 
neural representation (see peak in Fig. 5 C at HM). 

The location of the virtual visual streak, i. e., in- 
creased values of det E, depends on the particular 
form of VI. Our 3D surface model of VI averages 
over the observed substantial variability [51-[TU]. A 
statement about the location of the virtual visual 
streak independent of this variability can be done 
when we describe the location of the virtual visual 
streak relative to the curvature of VI. For this pur- 
pose one needs to consider the Gaussian curvature, 
because it is determined only by intrinsic properties 
of the surface. The Gaussian curvature of a point on 
a surface is the product of the maximal and minimal 
curvature of all curves passing through this point. 
As a consequence, a saddle-shaped surface has neg- 
ative Gaussian curvature while a spherical surface 
is positively curved. The entrance of CS is saddle- 
shaped and has therefore negative Gaussian curva- 
ture. The location where the fundus of the calcarine 
sulcus reaches its maximal depth has positive Gaus- 
sian curvature. The peak of the virtual visual streak 
is located at the transition between negative and 
positive Gaussian-curved (Fig. 5 D). While 

this provides only a rough estimate of the location 
of a virtual visual streak it is not only independent 
from the variability but can also be used as a pre- 
diction for selective magnification in other areas in 
the cerebral cortex. 



4 Discussion 

Cortical magnification is the local description of 
retinotopy, or in fact, of receptotopic maps in gen- 
eral. Tonotopic and somatotopic maps are other 
prominent examples [SHT] for which cortical magni- 
fication can be defined in the same manner as for 
the visual modality. Cortical magnification is the 
most fundamental quantity of cortical feature maps 
that describes how much cortex is devoted to pro- 
cess sensory information. A thorough mathematical 
foundation of magnification is needed, if we want 
to describe anisotropics and other symmetry break- 
ing constraints in such maps that may serve distinct 
purposes in sensory information processing. We will 



start our discussion with experimental evidence sup- 
porting our predictions, continue with methods able 
to verify them, and provide the context to current 
alternative tensor methods. Then, the limitation in 
symmetry properties of flat models of retinotopy and 
our assumptions concerning the 3D surface model of 
VI are discussed. We flnish with conclusions that 
can be drawn from our results on sulcal develop- 
ment. 



4.1 Experimental studies supporting meridional 
asymmetry 

The approximate form of human retinotopy has 
long been established from various methodologies. 
Some methods, which were mainly used before non- 
invasive imaging became available, obtain M indi- 
rectly by assuming a M-scaled quantity and mea- 
sured this. If A/ is known at various loci and also 
additional constraints are given, such as certain sym- 
metries, one can determine the layout of the global 
retinotopic map. Psychophysical studies that obtain 
an M-scaled quantity provide an elegant measure of 
cortical magnification mainly because they do not 
require an explicit reconstruction of the cortical sur- 
face. They greatly simplify the analysis bypassing 
any potential problems concerned with surface re- 
construction. Despite this advantage, these meth- 
ods suffer from limited resolution, poor quantifica- 
tion and are unable to directly measure M. They 
rely on the assumption of Af-scaling, a principle that 
could be violated. 

Anyway, evidence for the predicted meridional 
asymmetry is found in a study investigating percep- 
tual filling-in |31] and in a report about the shape 
of travelling migraine scotoma j32| . Both studies in- 
dicate that cortical magnification is not meridional 
asymmetric but selectively increased on or close to 
horizontal meridian (HM). 

The perceptual filling-in refers to the tendency 
of stabilized retinal stimuli to fade and become re- 
placed by their background. The time required for 
this illusion was investigated |31j . It varies with 
the azimuthal angle (/) of the target in the visual 
field. The filling-in was facilitated as the target po- 
sition changed from the horizontal to the vertical 
meridian. The time was maximal in the horizontal 
direction and minimal in the vertical. The abso- 
lute value differed in the several experiments, which 
were performed to exclude artifacts. But in general 
for a target presented at 8° eccentricity the response 
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Figure 5: (A) Arcal cortical magnification factor dot Af projected onto tire visual hemifield. The color bar 
indicates the normalized value in units of a? . The yellow blob at an eccentricity = 3° along the horizon 
indicates the virtual visual streak. (B) Linear cortical magnification M vs eccentricity 9 in units of a along 
the horizontal meridian at = 0° (HM) and the meridian at = 60° (dashed and solid lines on flat and 
curved model of VI, respectively). (C) Normalized surface area gain det E and (D) Gaussian curvature along 
HM and = 60°. 
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time for the horizontal condition was more than 30% 
longer than the one for the vertical. The author was 
aware of the fact that his findings can be explained 
by anisotropic cortical magnification. If magnifica- 
tion on HM is selectively increased, the neural repre- 
sentation of stimuli located on HM is increased too, 
and the time required for perceptual filling-in is pro- 
longed. 

Similar support for selectively increased cortical 
magnification on HM comes from the shape of trav- 
elling scotomta observed during a migraine attack 
with aura (Fig. 6 A). One of the first scientific reports 
on migraine with aura was given by Lashley |32j . He 
described in astonishing detail the spatio-temporal 
patterns of travelling scotoma in his visual field dur- 
ing migraine. The cortical magnification was not 
considered directly. But the shape of the scotoma in 
his left visual hemifield clearly hints to meridional 
asymmetry in M . Figure 6 A reproduces the spatio- 
temporal patterns drawn by Lashley. For the pur- 
poses of illustration it was modified: a polar grid was 
added, the first marked position of the scotoma is not 
shown, and the latest observed loci of the scotoma is 
shown as a curve being filled black. This particular 
shape and location was observed about 13 minutes 
after the onset of the attack. Two possibly under- 
lying neural patterns are shown in Figure 6 B and 
C. Assuming a constant speed of 2.5 mm/min of the 
pathophysiological process within VI, the scotoma 
at minute 13 (filled) is about 32.5 mm away from the 
neural representation of the fovea. Mapped with Eq. 
[3] and the parameter values a = 0.117 and b = 0.067 
this would correspond to 13.6° degrees eccentricity. 
This is in sufficient agreement with the position of 
the retinal blind spot marked by the dashed circular 
curve. 

Independent of the parameter values for a and b 
of the monopole map in Eq. Q , the scotoma shown 
in Figure 6 A cannot be mapped such that its gross 
shape changes. In other words, the narrow notch 
along the green path will not disappear. Therefore, 
either the width of the underlying cortical wave (a 
spreading depression [451146] ) varies along its path in 
the cortex (Fig. 6 B) or, much more likely, retinotopy 
is distorted from that of the monopole map (Fig. 6 
C). The latter is in agreement with our predicted 
selective increase of cortical magnification along the 
visual location of the steep fundus of the calcarine 
sulcus (green path) . This provides a test for our 3D 
model retinotopy. When the selectively increased 
magnification follows a path not exactly on HM but 



is shifted to the upper (or lower) field, the neural 
representation of this path should have the largest 
length from the neural representation of the fovea 
to most anterior part of the retinotopic map. If 
this path mainly overlaps with the steep fundus, the 
cortical representation of the horizontal meridian is 
shifted to the ventral (or dorsal) bank of calcarine 
sulcus. A qualitative analysis shows that one can 
measure linear cortical magnification precisely by 
comparable but simpler drawings of visual migraine 
aura [IJ . Unfortunately this data was pooled for all 
meridional directions. 



4.2 Alternative tensor methods 

By definition, methods that obtain M indirectly 
suffer from the possibility that the assumption of 
Af -scaling locally fails. Only direct methods that 
precisely delineate the retinotopic organization of 
human visual can measure cortical magnification. 
These technique are mainly used today to model 
structural and functional variation of the human 
brain. They can detect group differences, for exam- 
ple, disease effects on cortical anatomy [Ulil^, but 
also M was investigated [15l[l7l[50] . Problems occur 
in these studies by smoothing and flattening experi- 
mental data, because these techniques are based on 
mathematical models of retinotopy that introduce 
certain symmetries (see next Section). A recent ap- 
proach to estimate M [18] avoids many of the these 
problems because it is completely data-driven. It 
preserves geodesic distances by directly computing 
AI on the folded cortical surface. They found that 
M is similar in the dorsal and ventral compartments 
of VI within each hemisphere. A selective analysis 
of M on HM is not included. 

These direct methods to obtain M also use a ten- 
sor approach to morphometry. They are mainly re- 
lated to the Laplace-Beltrami operator. One may 
ask: How do these approaches compare to our ten- 
sor description? 

In these studies the speciflc cortical surface pa- 
rameterization is not essential, though much ef- 
fort has gone in the development of computation- 
ally advantageous parameterizations. If gyi denotes 
the cortical metric in an arbitrary parameterization, 
its components can be used to define the Laplace- 
Beltrami operator as a generalization of the Lapla- 
cian differential operator. In contrast, to derive the 
cortical magnification tensor M the parameteriza- 
tion is essential. The cortical surface parameteri- 
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Figure 6: (A) Four subsequent snapshots of a travelling migraine scotoma drawn by Lashley [32]. The smah 
cross at the left indicates the center of gaze, and the numbers state the time in minutes elapsed before the 
location of the scotoma was outlined. The point of time when the neurological symptoms were first noticed 
was chosen as the start. The scotoma travels across the left visual hemifield. (B) and (C) Two different 
neural patterns in VI that can explain the scotoma in (A). Either the noticeable notch in the scotoma pat- 
tern along the green path in (A) is reflected in the neural pattern by a narrow wave width along the neural 
representation of the green path (B) or the retinotopic map is distorted and the wave width is constant (C). 
In this case, the green path marks the fundus of the calcarinc sulcus. 
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zation by sensory coordinates is, in fact, the very 
idea behind cortical magnification. That is why we 
have introduced the term cortical expansion tensor 
C for this particular form of Qyi- For other ten- 
sor approaches it is natural and often convenient 
for the purposes of visualization to parameterize the 
cortex after inflating it to a sphere using the in- 
duced spherical coordinates. One has to define sul- 
cal curves as anchors, e.g., the Talairach coordinate 
system [51], to achieve that corresponding anatom- 
ical landmarks, such as the very cortical sulci, oc- 
cur at the same spherical coordinate across subjects. 
The Lagrangian formulation of a receptotopic map 
is a tensor approach that provides a natural and in- 
dependent parameterization of the cortical surface. 
According to this description data is pooled across 
subjects in the sensory reference frame, which sort 
of replaces anatomical anchors by functional labels. 
A similar idea is used by individual- subject region of 
interest analysis to identify functional labels . 

4.3 Influence of ocular dominance 

The wedge-dipole model of retinotopy |55] intro- 
duces two extensions two the monopole map. It is 
firstly introduced to model multiple fields (V1-V2- 
V3 complex) and also to improve the representation 
of peripheral data. The term "dipoZe" in its name 
refers to the latter, which is, however, for our study 
rather irrelevant. The corresponding mathematical 
changes that also introduce additional parameters 
are therefore not further considered here. We will 
only discuss the other extension referred to by the 
term "wedt^e" . In mathematical terms wedge refers 
to a methods that compresses the visual hemifield 
along iso-eccentricity lines (6'=constant) to a wedge 
with apex angle below tt. Usually, it is a half disc 
(apex angle equals 7r/2). This extension can actu- 
ally serve two purposes. Firstly, but as mentioned 
not relevant in the context of this study, the com- 
pression allows a unified model for the topography 
of the full visual field in areas VI, V2, and V3. Sec- 
ondly, the compression induces a topographic shear. 
Such a shear was stated to relate to the combination 
of two full representations of the visual field into 
VI by ocular dominance columns [56l[57]. The ori- 
entation of the stripes, in which ocular dominance 
organize and which define the direction of the to- 
pographic shear, may well relate to the 3D form of 
VI. 

In primates the patterns of segregated ocular 



dominance stripes run mostly parallel to the neu- 
ral representations of 9 coordinate lines (for the pat- 
tern in human see |58]). At least for no n- human 
primates in a region close to the fovea stripes tend 
to run horizontally |59| . In the direction perpendic- 
ular to the stripes, twice as much cortical space is 
needed to inject the two ocular representations of 
the visual field. Therefore, M depends on direction. 
For the sake of simplicity, consider the patterns of 
ocular dominance stripes runing parallel to 9 coor- 
dinate lines in the whole visual field. M is usually 
measured along the 9 direction perpendicular to the 
ocular dominance stripes. The value of AI should 
be only half as large in opposite direction parallel 
to the ocular dominance stripes, because only one 
visual hemifield is mapped along this direction. The 
effect of such a constant rotational shear, which is 
not accounted for in our model, is just a scale fac- 
tor of i. For example the fiat and curved model 
of VI would be half as wide. Since this induces a 
quasi-conformal map on the fiat surface model the 
symmetry arguments given above are still valid ex- 
cept for the scaling. However, in the foveal region 
stripes tend to run horizontally. The opposite scal- 
ing effect is obtained when M is measured only in 
this region. If data for M is obtained as an average 
from both regions the effects can partly cancel. A 
detailed analysis of the effect of dominance stripes 
on cortical magnification and cortical curvature is on 
open problem. 

4.4 3D model of VI and its parameterization 

We have made four partly interlinked hypotheses to 
construct the averaged 3D form of VI, namely that 
(i) two-thirds of the surface area lie within the CS 
walls, (ii) the ventral and dorsal half of VI are sym- 
metric with respect to the fundus, in particular, they 
extend equally long into the anterior direction, (iii) 
the shape of the walls can be approximated with 
smooth profiles (Gaussian-shaped profiles with con- 
stant standard deviation) and (iv) the shape of the 
fundus can also be approximated with smooth a half 
Gaussian-shaped profile. To obtain a curved retino- 
topic map on this surface area of VI, three further 
hypotheses are made: namely, that (v) the horizon- 
tal meridian maps to fundus of the calcarine sulcus, 
and (vi) that meridians are equidistantly spaced in 
the curved VI (Eq.[T2|), and (vii) in all other respects 
the retinotopy is goverend by the original monopole 
map. 
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The hypotheses (i)-(iv) are, to our mind, justi- 
fied because some smaller deformations within the 
gross shape of the calcarine suleus are probably not 
systematic and we assume that, on average, their 
influence is small compared to the curvature effects 
we investigate in this study. The course of calcarine 
sulcus varies widely among individuals and the lu- 
nate sulcus, which is roughly oriented vertically to 
it, appears to mark the anterior boundary of VI in 
non-human primates, but it is often missing in hu- 
mans [B^. Autopsy data suggest that VI proceeds 
farther anteriorly in the lingual gyrus [Sllinjj while 
fMRI data show that the dorsal and ventral com- 
partments of VI are at least similar in absolute ex- 
tent [18]. If VI proceeds farther anteriorly in the 
lingual gyrus this can be balanced by higher curva- 
ture of this gyrus. 

Our main results are robust to some small devi- 
ations of these hypotheses, because we only inves- 
tigate the influence of the gross sulcal pattern on a 
standard retinotopic layout. For example, we have 
replaced the Gaussian-shaped profiles (iii,iv) with 
cosine-shaped profiles with similar results. Further- 
more, if (v) holds only approximately and a meridian 
in the lower or upper quadrant of the visual hemi- 
field with a mild negative or positive shift away from 
horizontal meridian, respectively, maps onto the fun- 
dus, cf. |61] . we can still draw our conclusions as to 
where is the virtual visual streak. In this case, the 
retinotopy and in particular its derivation M can 
still be rather similar in the dorsal and ventral com- 
partments of VI, because hypotheses (v) and (ii) are 
interlinked and the ventral and dorsal half of VI (ii) 
can shift accordingly to respect the major symmetry 
between upper and lower visual quarter fields |18| . 

In fact, it provides a test to verify the predicted 
correlation between anatomical landmarks and func- 
tional properties if (v) does not hold exactly and in- 
stead of the horizontal meridian another path in the 
visual field, possibly with a variable azimuthal ofl[set 
with respect to the horizon, maps onto the fundus. 
Suppose, for example, that a meridional path with 
a positive shift away from the horizon (upper visual 
quadrant) is found with an increased cortical mag- 
nification, as suggested by data shown in Figure 6. 
Then the neural representation of horizontal merid- 
ian should be on the ventral bank, because accord- 
ing to our prediction this path is the location of the 
fundus in the visual field (if reversed retinotopically 
mapped). Overall, these hypotheses mainly estab- 
lish a symmetry between the neural representation 



of the upper and lower visual field quadrant, as found 
ill M- 

While the 3D form of VI affects the intrinsic 
properties of the surface, the layout of the retinal 
grid on this curved surface ultimately determines 
M. It was been chosen in a somewhat teleologi- 
cal manner by enforcing that meridians are equidis- 
tantly spaced (vi) while at the same time providing 
more cortex along the fundus that is created by the 
fold. Although this seems reasonable, one should 
notice that the location of the virtual visual streak 
depends on this assumption. Furthermore, because 
of this assumption, M is not sensitive to the precise 
shape of the profiles of the lift function (in particular 
the Gaussian-shaped profiles with constant standard 
deviation (iii)). An alternative and completely self- 
organized method can be provided by neural net- 
works such as a Kohonen net. To model retinotopy 
on the curved VI, the 2D net of the Kohonen layer 
must represent in its lateral connections the metric 
of the curved VI. The natural learning set is the 
cell density in the retina. We have investigated such 
networks and found similar results |37| . The major 
difference is that the increased areal magnification 
on the horizontal meridian is shifted slightly into 
the periphery, where the 3D form of VI has positive 
Gaussian curvature. Whereas the maximum areal 
magnification in this analytic study is found in be- 
tween the change from negative to positive Gaussian 
curvature (Fig. 5 D). 

4.5 Development of cerebral sulci 

On functional grounds, an increased cortical repre- 
sentation of the horizon with respect to the average 
cortical magnification on meridians (Fig. 5 A) serves 
a similar purpose as a visual streak. It is impor- 
tant to note that our model predicts where a vir- 
tual visual streak is located, while the existence of a 
virtual visual streak is a mathematical consequence 
of the fact that a surface with non-constant Gaus- 
sian curvature cannot be isometrically embedded in 
one with constant Gaussian curvature (e.g. intrinsi- 
cally fiat surface). Hence, the purpose of our tensor 
approach is to predict the shape and location of a 
virtual visual streak within VI. These are not in- 
trinsic properties of the cortical surface alone but 
also depend on retinotopy. The virtual visual streak 
therefore provides a link between cerebral sulci and 
encephalization of function because it links anatomy 
to retinotopy. The virtual visual streak is defined 
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for the visual modality as a selective increase of the 
number of cortical neurons mapped into the visual 
domain. It can in principal be defined more gener- 
ally for all senses and also for motor function. Such 
a quantity requires, as cortical magnification does, a 
notion of distance in the sensory or motoric space. 
For the sensory modalities like visual, auditory, and 
tactile, the distance between stimuli can be defined 
based on retina, basilar membrane, and body sur- 
face, respectively. The receptor densities must also 
be considered {dctE vs. detAf). For the sensory 
modalities like olfactory or also for kinetic modali- 
ties this concept is not readily clear. For example, 
given two odorants, is it possible to define a notion of 
distance between them? There is a stereotyped map 
of odorant receptor inputs in the olfactory bulb |62| 
but the concept of magnification requires two such 
maps. Nevertheless, we believe our approach can be 
applied to all modalities. It may even guide us in the 
search for a notion of sensory distance for modalities 
where it is not naturally given. 

Animals with a real visual streak have lower 
encephalization than human, even with rather 
lissenccphalic cerebral hemispheres p4H30] . When 
a retina has developed a visual streak, it is not pos- 
sible to map visual input spatially homogeneous onto 
a functional field of a lissenccphalic cortex without 
changing the global layout of the map. In other 
words, the M-scaling principle is likely to be locally 
violated. In such (e.g. macropus eugenii), 

a correlation between retinal ganglion cell density 
and M was found only along VM but not HM [27) . 
Hence, neural input is increased at the cortical site 
of the representation of the visual streak. A surface 
undulation can provide more cortical surface area for 
this increased input and thereby reestablishing the 
M-scaling principle. Dasyprocta leporina, another 
mammalian, in which M along HM approximately 
corresponds to its visual streak, shows gradual en- 
cephalization of the visual streak. Its retinotopy has 
nevertheless asymmetries that are not directly re- 
lated to the topography of the retinal ganglion cell 
density [25] . 

Irrespective of phylogenetic relationships, species 
that inhabit open spaces may independently have 
evolved a visual streak or a virtual visual streak 
that enables them to better scan the horizon. In 
other words, if we approach human retinotopy from 
the comparative standpoint, we can explain why the 
sulcal position is where it is, though not necessarily 
how sulci develop. The specific location of CS might 



be a convergent evolutionary solution to enhance vi- 
sual function on HM. The human retina poses a 
very mild visual streak |19II21| and even without tak- 
ing the visual streak into account, the A/-scaling is 
violated in all current theoretical models of human 
retinotopic maps, because they have a decreased cor- 
tical representation of the horizon with respect to the 
average cortical magnification on meridians. At the 
location where M-scaling is violated either less or 
more cortical surface area is available to process in- 
formation from a certain number of retinal ganglion 
cells. Therefore, our predicted selective increase of 
cortical magnification due to cortical folding, i.e., the 
virtual visual streak with selectively increased det E, 
indicates a form of encephalized visual magnifica- 
tion. If the virtual visual streak is actually a phylo- 
genetic shadow of a visual streak in human ancestors, 
which is supported by its mild occurrence [19 [ l21 j . 
thus it is not an independent evolutionary solution, 
a switch from selectively increased retinal to cortical 
magnification on HM took place. 

To summarize, the theme from structure to func- 
tion is central to biology, as is the reverse direction 
from function to structure. We doubt that the pre- 
cise correlation between the primary visual cortex, 
as an functional area, and the calcarine sulcus, as an 
anatomical structure, occurred by chance. 

None of the existing mathematical models is 
suited to describe retinotopy in relation to the cor- 
tical geometry. Most of them are in fact either ID 
functions assuming invers linear mapping given by 
two parameters {a,b or A,E2, see Eq.[3]), without tak- 
ing into account that a 2D generalization is not triv- 
ial, or they are intrinsically flat maps described by an 
analytical function. The monopole map and its five 
parameter extension, called wedge-dipole map [55] 
are by definition conformal and quasi-conformal, re- 
spectively. Any analytic function has a high inner 
structure, i.e., exhibit certain symmetries. Alterna- 
tive nonconformal maps suffer from inconsistencies 
with the definition M [G^ . It is assumed that single 
components of the Jacobian matrix J$ (Eq. [5|) scale 
with M. But J$, unlike M, contains a rotation 
carrying the retinal directions onto its neural repre- 
sentations. This renders the assumption implausible 
if conformal mapping is not presumed in the first 
place. 

Generally, when cortical functions became en- 
cephalized and even completely new sensory and mo- 
toric areas simultaneously arose, this constitutes a 
developmental force for the prominent patterns in 
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the highly convohitcd surface of the cerebral cor- 
tex in human. Linking cortical magnification to 
anatomy with methods of continuum mechanics and 
complex analysis may shed some light on what could 
have been the underlying pattern formation princi- 
ples. 
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